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Abstract 

A triple vector correlation in the /i+ e+e+e~ decay with polarized 
muons is investigated as a probe to CP violating coupling constants 
in supersymmetric models. A sizable triple correlation can be induced 
due to a complex phase in the supersymmetric soft-breaking terms in 
the SU(5) grand unified theory. Correlation with the electric dipole 
moments of electron and neutron are investigated and it is shown that 
these quantities give independent information on possible CP violating 
sources. 
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Lepton flavor violating (LFV) processes such as /i+ —>■ e+7, /i+ e^e^e~ and 
fi~-e~ conversion in atoms have important implications in search for physics beyond 
the Standard Model (SM). Many extensions of the SM predict measurable rates 
for these LFV processes. In particular it has been pointed out that rates for these 
processes can be as large as or just below the present experimental upper bounds in 
supersymmetric (SUSY) grand unified theories (GUTs) Current experimental 
bounds for these processes are B(yU+ < 4.9 x lO"^"*^ [Q, B(/i+ e+e+e^) < 

1.0 X 10-12 a{fx-Ti e-Ti)/or(/i-Ti capture) < 4.3 x lO^^^ g_ Further 
improvements of these bounds in two or three orders of magnitude will have a great 
impact on search for unified theories based on SUSY. 

In this letter, we consider the /i+ e^e~''e~ process with polarized muons. If 
the initial muon is polarized we can define T-odd triple vector correlation, a- {P1XP2) 
where a is the muon spin and pi, p2 are two independent momenta of the final 
particles. In this way we can measure CP and T nonconserving effects in LFV 
interactions ||^. As possible sources of CP violation we consider SUSY soft-breaking 
terms in SUSY GUT. We show that it is possible to generate a measurable T-odd 
asymmetry in this model. We also consider the electron and the neutron electric 
dipole moments (EDMs) which are sensitive to the same CP violating phase and 
show these observables give independent information on the possible CP violating 
sources. 

Let us begin with the following effective lagrangian relevant for the /^^ —>■ 
e~^e~^e~ decay. 

+ 9i eRCL + 92 W^^R eZe/j + g^i JlR'y^eR e^'jf.eR 

+ 94 Wl'^eL eZ'j^eL + 95 Mi'^^r cZi^cl + 96 Wi'^^l e^T^eiJ + H.C. (1) 

Here Ar, Al and Qi^s are in general complex coupling constants. In order to obtain 
the above expression we have used the Fiertz rearrangement for four-fermion terms 
and neglect terms suppressed by the electron mass compared to the muon mass. 
Among the above terms the photon-penguin terms, Ar and A^, induce the 
—* e^7 decay in addition to the /i"*" — >■ e~^e~^e~ decay. The branching ratio 
for /i"*" — i> e"'"7 is given by B(/i"'" e~^"y) = ^i- + On the other 
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hand the /i — > 3e branching ratio depends on four-fermion couphng constants gi^s 
as well as and Al. In order to present the differential branching ratio for this 
process we first discuss fi ^ 3e kinematics with polarized muons which is specified 
by two energy valuables and two angle variables. We assign p, pi, p2 and ps to the 
momenta of fi~^ and two e+'s and e~, then two energy variables are taken to be 
Xi = and X2 = ■ Two angle variables are necessary to specify the relative 
position between the muon polarization direction {n) and the decay plane. Defining 
Pi for the momentum of the positron which satisfies (n x ^3) -pi > 0, angle variables 
{6, (f)) can be identified as a polar coordinate of n (0 < 6* < tt, < < tt) in 
the coordinate system where is taken as the z direction and the decay plane is 
identified as z-x plane with {pi)x > (Fig. |l]). Note that in this notation the T-odd 
asymmetry appears as an asymmetry in exchange of Xi and X2 in the Dalitz plot 
after integrating out the angle variables. 

It is now straightforward to calculate the differential branching ratio from the 
effective lagrangian. A detail formula will be given elsewhere [^. There are three 
types of terms in the differential decay branching ratio, i.e. square of photon-penguin 
terms, square of four-fermion terms and interference between the photon-penguin 
terms and four-fermion terms. T-odd terms arise as a part of the interference terms. 
These are given by 



X [ 2(xi +X2- l)Im(^3^I + g^A*^) - (2 - xi - X2)lvci{g^Al + ^e^n) ] ,(2) 



where P is the polarization of muons and e is the positron charge. As seen from 
the above expression the T-odd term is generated if the photon-penguin term A^ 
(or Afj) has a different phase from the four-fermion terms g^, g^ (or (74, g^). These 
terms are in fact odd in exchange of Xi and X2. 

In order to extract these T-odd terms from experiments we also need to know 
distribution of other terms in the /i ^ 3e differential branching ratio. In particular, 
the \Ai\'^ and terms have the following form. 

3 

c^Bphoton = ^2 dxidx2dcos 9d(j) 



3 



(ia;i(ix2(i cos 6'(i0P sin 6* sin 8e(xi — X2) 




2 



8e' 



2x\ - 2xi + 1 ^ 2x\ -2x2 + 1 



(|A^j2 + |Az.r), (3) 



1 — X2 \ — X\ 

where we only keep terms independent of the muon polarization and neglect the 
electron mass. The above expression is singular for x\^\ and X2 — * 1 if we neglect 
the electron mass, and therefore the integrated branching ratio is solely dominated 
by this kinematical region if eAi^ and eA^ are similar in magnitude as (/j's. In such 
a case the relation B(yU — > 3e)/B(/i e'j) ~ 1/150 is known to hold 0. Since our 
purpose is to look for the interference between the photon-penguin and the four- 
fermion terms we should exclude the region very close to Xi = 1 and X2 = 1. In the 
actual experiment T-odd asymmetry should be extracted from the investigation of 
distribution in the Dalitz plot. Here in order to present the T-odd effect we use the 
branching ratio and asymmetry integrated in the following way. Defining the regions 
Ri = {xi+X2 > 1, X2 < X2 > Xi} and R2 = {xi+X2 > 1, Xi < l—S, Xi > X2}, 
where 6 is introduced to cut off the region Xi,X2 — 1, the integrated branching ratio 
and the T-odd asymmetry are defined as 

B[S] = f dxidx2 C dcose r d /^}^] ^ (4) 



R1+R2 J-i Jo dxidx2dcost 

^ J n r ^ e+e+e^) 



/i r7i 
dcosO J 



PB[6] 

dxidx2 j dcosO J 



Ri J -I Jo dxidx2dcos 

1 ,rfB(/i+ ^ e+e+e^ 



^^^^^^ /-www n.j. ■ (5) 

IR2 J -I Jo dxidx2dcos 

Let us consider the SU(5) SUSY GUT model. In this case the LFV occurs 
from the loop effect between the Planck scale and the GUT scale. Even if we assume 
that all scalar fields have a common SUSY breaking mass at the Planck scale, the 
large top Yukawa coupling constant becomes a source of the flavor mixing in the 
slepton sector since sleptons belong to the same GUT multiplet as squarks |^, ^ . 
In the simplest SU(5) SUSY GUT the Yukawa couplings are given by the following 
superpotential. 

W = T,{fuhT,H + T,{fdh,F,H, (6) 

where T, are 10 dimensional representations and Fi are the 5 dimensional represen- 
tations and H{H) is 5(5) dimensional Higgs superfield. Due to the loop effect of 
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fu the right-handed stau becomes hghter than other right-handed selectrons at the 
GUT scale and therefore slepton mass matrix is no longer simultaneously diagonal- 
ized with lepton mass matrix. In the approximation that the first two generation's 
sleptons are degenerate every term in the LFV amplitude is proportional to Xr = 
^R(e)Te^i?(e)r/x where Vji{e) is a right-handed lepton mixing matrix in the bases 
where the slepton mass matrix is diagonalized up to the left-right mixing terms 
1^. Therefore if there are no other source of complex coupling constants almost no 
asymmetry is generated because the photon-penguin and four-fermion terms have 
approximately the same phase. Situation will change if we allow complex phases 
for the SUSY breaking terms. Within the assumption of the universal soft-SUSY 
breaking terms we can introduce two independent phases which we take a phase of 
the trilinear coupling constant {A term) and a phase of the higgsino mass term (/i 
term). These phases in general induce large EDMs of neutron and electron if the 
masses of squarks or sleptons are in the range of a few hundred GeV |T^. In this 
respect an interesting observation was done in Ref. |TT| that unlike the fi phase, the 
constraint on the A phase is much weaker so that even a 0(1) phase is allowed. In 
the foUowings we assume that the fi phase vanishes for simplicity and see effects of 
the A phase to the electron and the neutron EDMs and the T-odd asymmetry in 
the /i — i> 3e decay. 

We calculate the branching ratio Eq.(H) and the T-odd asymmetry Eq.(^ 
and the electron and the neutron EDMs in the SUSY SU(5) GUT model with a 
complex A parameter. For calculations of the electron and the neutron EDMs in 



SUSY models we follow Ref. . Assuming the universal soft-breaking terms at the 
Planck scale we solve the renormalization group equations for the coupling constants 
and SUSY soft-breaking terms from the Planck to the GUT scale, and then the GUT 
to the weak scale. We also require that the electroweak symmetry is broken properly 
due to the renormalization effect and take into account phenomenological constraints 
from various SUSY particle searches including the constraint from the 6 — > 57 decay 
as discussed in Ref.|]13[. Once we obtain the SUSY breaking terms at the weak 
scale we can evaluate the coupling constants A^, Al and (^j's in Eq.([l|) through 
one loop SUSY diagrams. Ar, Al are determined by photon-penguin diagrams and 
QiS get contributions from off-shell photon-penguin, Z-penguin and box diagrams. 



Formulas for these contributions can be found in Ref. 
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In the calculation of the /i — > 67 and — > 3e branching ratios there is an 
important ambiguity associated with the Yukawa coupling constants at the GUT 
scale. As discussed before the branching ratios are proportional to |Arp. If the 
Yukawa coupling constants are given only by Eq.(^ the lepton mixing matrix can 
be related to the quark's Cabibbo-Kobayashi-Maskawa (CKM) matrix element. It 
is known, however, that this assumption does not lead to a realistic mass spectrum 
for fermions. As discussed in Ref. if other operators are relevant for generation 
of Yukawa coupling constants below the GUT scale the relationship between the 
lepton mixing matrix and the CKM matrix becomes quite model-dependent and A,- 
can be significantly different from the value obtained with the above assumption. 
In view of this ambiguity we treat as a free parameter and calculate the /x — >■ 67 
and fi —>■ 3e branching ratios normalized by lA^-p. Note that the asjTiimetry and 
the EDM are essentially independent of A,-. 

In Fig.|^. (a) we show the branching ratio for /i — 67 and /i — 3e normalized by 
I At- p as a function of the right-handed selectron mass [m-^). The SUSY parameters 
in this model are taken as the SU(2) gaugino mass M2, a complex A parameter 
defined as Csoft = —maAx{Ti{fu)ijTjH + Ti{fd)ijFjH) at the Planck scale, and 
the ratio of two Higgs vacuum expectation values (tan (3) and the universal scalar 
mass mo at the Planck scale and the sign of fi. In this figure we fix M2 = 200 
GeV, Ax = i and tan/5 = 3,10,30, fi > and the branching ratio are shown 
as a function of instead of mo. The cut off parameter 6 is taken to be 0.02 
for the /i — s> 3e branching ratio. We fix the top quark mass as 175 GeV. If, for 
example, |At-| = 1 x 10^^ the /i — > 67 branching ratio is 10~^°-10~^^ and the /i ^ 3e 
branching ratio is lO^^^-lO^^^ for tan f3 = 10. On the other hand if At is given by the 
corresponding CKM matrix elements Vi^ts then |At-| = (3 — 5) x 10"'^ and therefore 
the branching ratio is smaller by three orders of magnitude. In Fig.|^.(b) the T-odd 
asymmetry defined in Eq.(|]) is shown for the same parameters as in Fig.|^.(a). Also 
the electron and the neutron EDMs are shown in Fig.0.(c). We can see that the 
asymmetry becomes maximal around m~^ = 400 GeV while the magnitude of the 
electron (neutron) EDM is a decreasing function of m~^ above 200 (400) GeV. The 
difference of these behaviors comes from the fact that for the asymmetry photon- 
penguin and four-fermion terms have similar magnitude for m~^ ~ 400 GeV but 
for the EDMs only photon-penguin diagram is relevant. It is interesting to see 
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that when the asymmetry is large the ratio of the B(yU 67) and B(yU — > 3e) 
is 1/50 - 1/100. Fig.|^.(a) shows the correlation between the absolute value of the 
electron EDM and the T-odd asymmetry. In this figure, the SUSY parameters are 
scanned in the region < mo < 2 TeV, \Ax\ < 5, < Mx < 2 TeV and fix 
aig{Ax) =vr/2, /i > 0. We find that the sign of the EDM and that of asymmetry 
are the same but here the absolute value of the EDM is shown. We can see that the 
electron EDM and the T-odd asymmetry do not show strong correlation. Within 
the current experimental upper bound of the EDM, which is \de\ < 4 x 10^^'' |jl5 



the asymmetry can be as large as 18 %. This means that the T-odd asymmetry in 
the n 3e process gives independent information from the electron EDM on the 
possible complex parameters in the SUSY GUT model. Correlation between the 
neutron EDM and T-odd asymmetry is also calculated. As in Fig.^(a) we do not 
see any correlation and the ±18% asymmetry is possible within the experimental 
upper bound of the neutron EDM (|c/„| < 1.1 x lO^^s Fig.|.(b) shows that 

the correlation between B(/i ^ 3e)/|AT-p and T-odd asymmetry. In this figure we 
take into account the experimental bounds of the electron and the neutron EDMs. 
We see that, for example, when |A| is 1 x 10^^ 10% asymmetry is possible for 
B(/x 3e) ^ 10"^^. We also investigated the correlation between the asymmetry 
and the slepton mass and found that the slepton mass which corresponds to maximal 
asymmetry changes from 200 GeV to 2 TeV when we scanned the parameters in the 
above region, so that a large asymmetry does not necessarily mean a light slepton. 

Let us discuss T-odd asymmetry in SUSY models other than SU(5) GUT. 
In the minimal SO(IO) SUSY GUT model it is pointed out in Ref. p|] that the 
/i — s> 67 branching ratio is enhanced by (mT-/m^)^ compared to the SU(5) model. 
The /X — s> 3e branching ratio is also enhanced by the same amount and in this 
case the photon-penguin and four-fermion terms can have different phases without 
the phases of the soft breaking terms. But unfortunately, the T-odd asymmetry 
cannot be large because only the photon-penguin term is enhanced in this decay. 
Another possibility to induce large LEV is the SUSY model with the right-handed 
neutrino supermultiplet |]14|, [1^. In this case the right-handed neutrino's Yukawa 
coupling constants become new sources of LEV and CP violation. This model is 
also interesting in connection with the baryon number asymmetry of the Universe 
since the right-handed neutrino or right-handed sneutrino decays can be the origin 
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of the leptogenesis |T^ and CP violation in the Yukawa couphng constants is one of 
the necessary ingredients. We have investigated the LFV branching ratio and the 
T-odd asymmetry in this modeL We take neutrino masses in the range of 10^^ - 10 
eV and adjust the right-handed neutrino Majorana mass to 10^^ - 10^^ GeV to get 
0(1) Yukawa couphng constants. In this model using the freedom of right-handed 
Yukawa coupling constants, it is possible to generate two independent phases in 
the photon-penguin terms and the four-fermion terms without help of the complex 
soft breaking terms. With real soft-breaking terms typical magnitude of asymmetry 
turns out to be less than 0.1 % although it is possible to have an asymmetry up to 
10 % by tuning the parameters in the right-handed Yukawa couplings. 

We only consider the polarized muon decay here. Extension to the tau decay 
is straightforward. For example the T-odd asymmetry in r ^ 3yU is obtained by 
replacing relevant generation indices in the fi —>■ 3e formula. It should be noticed, 
however, that there is essentially no difference between SU(5) and SO(IO) models 
in this case because there are no enhancement mechanism for the SO (10) model 
compared to the SU(5) case in this decay mode. 

In this letter we investigate a possibility to observe a sizable CP nonconserving 
asymmetry in the fi —>■ 3e decay in SUSY models. In the SU(5) GUT model sizable 
LFV arises from the GUT interactions and an asymmetry up to 18 % can be induced 
by the complex phase of trilinear soft-breaking term. In order to measure the asym- 
metry of this magnitude the the /x ^ 3e branching ratio has to be large enough. 
In the SU(5) model the magnitude of the branching ratio itself strongly depends 
on the model-dependent parameter so that practically the fi 3e branching 
ratio is only constrained by the experimental bounds of the /i — > 3e and ^ ^ e'y 
processes. Therefore if the /x — 67 process is observed at the level of lO^^^-lO^^^ 
in near future, the /i — > 3e branching ratio can be 0(10^^'^) and an experiment of 
fi 3e with a sensitivity of the 10~^^ level could reveal the T-odd asymmetry. Note 
that stronger bounds on EDM or slepton mass do not mean that the asymmetry 
cannot be measurable because due to ambiguity of A,- the branching ratio cannot be 
strongly constrained by these bounds. Besides the EDMs for neutron and electron, 
the T-odd asymmetry therefore gives us a new possibility to search for CP violation 
in SUSY models. 
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Figure Captions: 



FIG. 1: Kinematics of — > 3e decay in the center-of-mass system of muon. The 
plane I represents the decay plane on which the momentum vectors pi, p2, ps lie, 
where pi and p2 are momenta of two e~'"'s and ps is momentum of respectively. 
The plane II is the plane which the muon polarization vector n and p^ make. 



FIG. 2: (a). Branching ratios for fj, ^ and // — > 3e normalized by |A^p= 
|l^R(e)Tel^(e)T/iP as a function of the right-handed selectron mass m~^. The cut-off 
parameter S is taken to be 0.02. We fix the SUSY parameters as M2 = 200 GeV, 
Ax = i, fJ^> and tan/3 = 3 (dotted line), 10 (solid line), 30 (dashed line) and top 
quark mass as 175 GeV. (b). T-odd asymmetry A as a function of m~^. Parameters 
are the same as in (a), (c). The electron and the neutron EDMs as a function of 



FIG. 3: (a). A correlation between T-odd asymmetry A and the absolute value 
of the electron EDM \de\- The SUSY parameters are scanned in the region < 
mo < 2 TeV, \Ax\ < 5, < Mx < 2 TeV and fix arg(Ax) = 7r/2, tan/3 = 5 and 
/J, > 0. (b). A correlation between T-odd asymmetry A and the branching ratio for 
^ 3e normalized by jArp. The SUSY parameters are the same as in (a) and the 
experimental bounds of the electron and the neutron EDMs are imposed. 
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Fig.l 
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Fig.2.(c) 
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Fig.3.(b) 
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